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1. Introduction 



The ternary semiconducting compounds A^^ Bl^^C"^^ have been widely investigated 
because of their potential applications in electro-optic, optoelectronic, and non-linear 
optical devices Most of these compounds have the defect chalcopyrite or stanite 
structure p|, 0]. In an A^^ CX^ defect chalcopyrite compound. A, B, and C atoms 
and the vacancy E, are distributed as follows 0,: A site at (a/2,0.c/4), B at (0,0, c/2) 
and (0,a/2,c/4), C at 7), (a, /5, 7), (/?, 0,7) and (/?,«, 7) and E at (0,0,0). It's 

lattice parameters are given by a = (3 = a/4,7 = c/8 and c = 2a, as shown in the 
figure 1. 

In this communication we shall investigate the electronic structure and optical 
properties of the defect chalcopyrite ZnIn2Te4 from a first principles approach. The 
chosen system is not only important for its application in technology, but is also 
important for the study of the effect of defect or impurity on various properties of 
materials. There is also growing interest in the study of defect and nano structures 
1^ , . One of the experimental method for the preparation of nano-structures is to fill 



in the structural voids in a "host" compound with atoms of a given substance. These 
structures could either be local void clusters or even tubulary voids. The characteristics 
dimensions. The host compound helps to form a matrix in which these nano structures 
stabilize. In a subsequent communication, we shall extend the present work to include 
nano structures in host materials. 

Recently Ozaki et al [^, ^ have carried out a detailed experimental and theoretical 
study of optical properties of amorphous and crystalline ZnIn2Te4. There are few other 
experimental measurement [§- [l^ on the compound also but there has been no first 



principles calculation for the electronic structure and optical properties of this material. 

The band structure calculations by empirical parameterized tight-binding methods 
has been carried out for ZnIn2Tei by Meloni et al [jlT| and Ozaki et al These 
calculations required fitted parameters. The actual crystal structures also seem to have 
been simplified. Meloni et al have assumed a pseudocubic structure (space group = 
Vd), rather than the actual chalcopyrite type structure. Ozaki et al have taken a defect 
chalcopyrite structure. However, the space group Si, which they assumed, does not 
reflect the correct symmetry of this structure. This is the usual chalcopyrite structure 
but with vacancies at the sites as shown in the figure 1. Because of the vacancies, this 
structure does not have the full space group Si- Rather, we can assign this structure the 
/4 symmetry. The positions of Te's are zinc-blende type, i.e. two inter-penetrating fee 
lattices shifted one-fourth of the way along a body diagonal. The values of the lattice 
parameters a and care taken from experiment to be 6. 11 A and 12.22A respectively as 
reported by Hahn and Ozaki |^. The positions of the atoms within the unit cell is 
shown in table 1. 
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Table 1. Positions within the unit ccU of atomic basis, including the empty spheres 
to take into account the voids within the structure 



2. Results and Discussion 

2.1. Electronic Structure Calculations 

To start with, for our electronic structure calculations we have used the well established 
tight-binding linearized muffin-tin orbitals method (TB-LMTO), discussed in detail 




Figure 1. Crystal structure of defect-chalcopyrite-type semiconductor ZnIn2Te4. 
Here an orthorhombic primitive cell is shown 



Electronic and Optical Properties of Znln2 Te^ 



4 




z r X p N 



Figure 2. The band structure of ZnIn2Te4 within the TB-LMTO 



elsewhere lib, 17, O. Electron correlations are taken within local density approximation 



of density functional theory [|1^, ^ 

The TB-LMTO energy bands in several high symmetry directions in reciprocal 
space is shown in figure Let us examine this in some detail. The basis of the TB- 
LMTO starts from the minimal set of muffin-tin orbitals of a KKR formalism and then 
linearizes it by expanding around a 'nodal' energy point E^^. The wave-function is then 
expanded in this basis : 
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CI and A2 are TB-LMTO potential parameters and S'££j(k) is the structure matrix. 

First note that for the calculation of the optical properties of the solid we need 
to span a large energy range, from the occupied valence to the unoccupied conduction 
states. For ZnlnTe this spans a range from -15 eV to 10 eV. The "nodal" energies 
cluster below -5 eV in the valence band. Obviously, the conduction part of the band is 
not accurately reproduced. The third generation NMTO expands the wavefunction in 
term of a basis which is expanded as a Lagrange interpolation around a discrete set of 
nodal energies {e„} : 



N 

(2) 



E CL(r) 4%(k) 



n=0 



The C^nRL R'L' Lagrange matrices which are such that the energy dependent partial 
wave basis (j)1{r,E) takes the values 4>ni{r) at the nodal energies. Unlike the LMTO, 
the nodal energies are independent of the indeces La. By choosing the nodal energies 
across the energy range of interest we can accurately reproduce the bands in that range. 
This was shown for GaAs in the range -15 eV to 20 eV by Andersen and Saha-Dasgupta 



The figure ^ shows the energy bands obtained from the NMTO using three nodal 
energies spread across the energy range of interest. In comparison with the figure 0, we 
note that the largest change occurs in the conduction band, away from the TB-LMTO 
nodal energies bunched below -5 eV. We have used the von Barth-Hedin exchange 
with 512 k-points in the irreducible part of the Brillouin zone. 

The calculated density of states within TB-LMTO and NMTO are shown in the 
figure It is evident from these two figures that there is significant difference between 
these two only in the conduction band. From the band structure comparison, this was 
also evident for the reasons discussed earlier. These densities agree quite well with the 
experimental XPS measurements of Ozaki et al and also with their parametrized 
tight-binding calculations. These results are shown in the figure |. It must be noted 
here that since these calculations are based on the LDA, we do obtain a lower band 
gap. In this case the NMTO estimate of the band gap is around 1.3 eV, while the 
XPS data indicate a gap of around 2 eV. Nor can we have much confidence about the 
conduction bands. This problem can be tackled more accurately, for example, through 
by quasi-particle band structure within a GW type approximation |2^. These NMTO 



based calculations are a much better starting point of GW self-consistency iterations 



than the first or second generation LMTOs (see comments in |TB]) 

The TB-LMTO band structure gives us insight into the various structures in the 
density of states. At the lowest energies around -14 eV we have the states arising out of 
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Figure 3. The band structure of ZnIn2Te4 within the NMTO 



the Te-s electrons. The sharp peak at -10 eV is due to the very narrow Zn-d bands. The 
next structures around -8 eV arise from the In-s and Zn-s electrons. The predominantly 
covalently bonded Te sp states gives rise to bonding and anti-bonding bands around -5 
eV and 5 eV respectively. The band gap lies between these bands. 



2.2. Optical Properties: 



In recent years a number of methods have been proposed for calculating optical 
properties within the framework of the LMTO ||2l|, ^ for both metals [0, ^ 

and semi-conductors ||2^ 



^ |, Upenski et al |^ proposed a method for the 

accurate calculation of the optical matrix elements based on the continuity equation 
for the charge -density operator. They were able to calculate accurate optical matrix 
elements by including the combined correction term to compensate for inaccuracies in 
the wavefunction due to the basis set in LMTO theory being finite. This method and 
corrections were subsequently applied by other authors ^ to calculate accurate 
optical spectra of Ge, GaAs, InSb and CdTe. 
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Figure 4. Density of States for ZnIn2Te4 (top) calculated from the TB-LMTO and 
(bottom) calculated from NMTO with three nodal energies across the energy range 



In all the above calculations however, the gradient operator has been employed for 
the determination of the optical matrix. There is an alternative method developed by 
Hobbs et al which avoids the determination of gradient operator. Their method 
allows for the inclusion of non-local potentials in the Hamiltonian. In their method they 
employed Green's second identity and the commutation relation between the position 
and Hamiltonian operators. They finally wrote the momentum matrices in terms of 



Gaunt coefficients |28[ and potential parameters which are defined within the LMTO 
method [|16[ The mathematics has been described in detail both by Hobbs et al 
and by Saha et al ([0). We shall indicate here the minor changes required for the 
calculation of the transition matrix element within the NMTO. 
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Figure 5. (a) Theoretical density of states N(E) and (b) experimental XPS spectrum 
1(E) taken from Ref. [| 



The expression for the imaginary part of the dielectric response remains the same, 
as derived from the Kubo formula : 

'^^"^^ = ^^^^T.T.me.-mfF.i'^-FMEf-E^-hu) (3) 

where, m* is the effective mass of the electron, f2 is the volume of the sample, i and f 
refer to the initial and final states respectively, 7 refers to the direction of polarization 
of the incoming photon, p is the momentum of the electron and Fij is the occupation 
probability of the initial and final states respectively. For semiconductors, i lies in the 
valence band and f in the conduction band and at O^K we have Fi = 1 and Ff = 0. 

We can obtain the real part of the dielectric function ei (u) from a Kramers-Kronig 
relationship. 

2 (uj' - u) €2(00') , , 

TT JO Ul' — Ul^ 
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The transition probability which involves momentum matrix elements can be 
evaluated using a gauge independent formalism and commutation relation 



dv ruer . 



Thus calculation of momentum matrix elements means calculation of the following 
integrals: 



By noting the relations: 

H rndr) = el rndr) 
and using Green's second identity we may write the integrals as: 

/ CL'(r)^ r H ^nd^) dh = e: T^l' fj €,{r) 4>l,{r) dr 
J CL'lr)-^ H r 0^Jr) dh = /-''T^l' £ r^'ir) Keirydr.. . 



where, Sq, is the atomic sphere radius of the a-th atom in the unit cell and is a 
combination of Gaunt coefficients : 

T,,, = [ (g-[q''- - Gtir) i + 1 + G%^f) j + v~2 Gtfr t 

D^i^ is the logarithmic derivative of at r = Sq, 

The lattice for ZnIn2Te4 does not have cubic symmetry. If we take the c-axis to 
be the z-axis, it is clear that the three optical responses e^, €y and are not the same. 
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Figure 6. (Top) Real and Imaginary parts of the dielectric function of ZnIn2Tei 
calculated from an NMTO calculation. The full lines show e±_ and the dashed lines e||. 
(Bottom) The experimental results of Ozaki et al 

In the top part of figure ^ we show the variation of the real and imaginary parts of e\\ 
= ez and e_L = e^^ + with the incident photon energy 

We first note that since our LDA-based calculations give a smaller band gap that 
experiment, we have applied the scissors operator, which involves a rigid shift of the 
conduction band with respect to the valence band, so that the band gap matches. This 
is frequently used by LDA practitioners, but cannot be fully justified. The correct 
procedure would be to carry out a full GW calculation which gives rise to an energy 
dependent self-energy, which shifts the bands unequally at different energies, resulting 
in a distortion of the shape of the densities of states as well. Given this, the agreement 
of the theoretical calculations with the experimental results of Ozaki et al ^ available 
up to 10 eV photon energies, is not bad. The experiment does not align the crystal 
axis with the polarization of the electric field, so that we should compare the results the 
direction averaged response (ey + e±)/3. The principal peak positions and heights are 
well reproduced. 
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We shall conclude with the remark on the two directions which we propose to 
take from here. The first point is to recognize that the NMTO calculations form a 
reasonable starting point of the more sophisticated many-body GW approaches. The 
second is related to the reason why we chose to study the defect chalcopyrite in the 
first place. We would like to fill the voids in the structure with various "impurities" 
and study the signal for these in the optical response. These will be the aim of our 
subsequent work in this area. 
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